Decoherence in a scalable adiabatic quantum computer 
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We consider the effects of decoherence on Landau-Zener crossings encountered in a large-scale 
adiabatic-quantum-computing setup. We analyze the dependence of the success probability, i.e. the 
probability for the system to end up in its new ground state, on the noise amplitude and correlation 
time. We determine the optimal sweep rate that is required to maximize the success probability. 
We then discuss the scaling of decoherence effects with increasing system size. We find that those 
effects can be important for large systems, even if they are small for each of the small building 
blocks. 



I. INTRODUCTION 

The promise of enormous levels of speed up over classi- 
cal computing algorithms has stimulated research in the 
field of quantum information processing, especially af- 
ter the discovery of a variety of concrete algorithms, in- 
cluding the factoring and search algorithms In the 
commonly studied approach, to which we shall refer as 
sequential quantum computing (SQC), the calculation is 
performed using a sequence of pre-designed unitary op- 
erations on the quantum state of the system. An alter- 
native to SQC was proposed a few years ago, namely 
adiabatic quantum computing (AQC) @, The main 
motivation for pursuing AQC is the idea that certain cal- 
culations could be performed with speeds comparable to 
those obtainable with SQC using a drastically different 
approach that avoids some of the difficulties associated 
with SQC. 

Calculations in AQC are performed as follows: one 
takes a given quantum system and sets the external pa- 
rameters such that the system is guaranteed to relax to 
its ground state. One then slowly varies those external 
parameters until the desired final set of parameters is 
reached. The result of the calculation is then encoded 
in the final quantum state, which should be the ground 
state of the Hamiltonian at the end of the process. Dur- 
ing this adiabatic variation of parameters, a large num- 
ber of avoided level crossings are encountered, and the 
physics of Landau-Zener (LZ) transitions applies 0, @. 
The LZ formula, which will be given below, states that 
if the time taken to sweep across an avoided crossing is 
long compared to the inverse of the gap in that crossing 
(we take h = 1), the system remains in its ground state 
with a high degree of certainty. 

The fact that in AQC the system remains in its ground 
state suggests, at least at first sight, that AQC is robust 
against decoherence H, 01- In fact, that robustness is 
generally thought of as being the single major advan- 
tage over SQC. Recently it has been argued, however, 
that decoherence does set limitations on AQC @, H, @ • 
In particular, if the passage from the initial to the final 
state is done too slowly, the success probability of the 



algorithm will be reduced from the maximum obtainable 
value. In this paper we analyze the optimal implemen- 
tation of an AQC algorithm in the presence of a noise 
source. We also discuss how decoherence effects increase 
in importance with increasing system size. We show that 
decoherence considerations can play a major role in de- 
termining the optimal operation conditions of a scalable 
AQC setup. 

This paper is organized as follows: In Sec. [Til we 
present the basic LZ problem. In Sec. IIIII we briefly 
comment on the question of the scaling of the minimum 
gap with system size. In Sec. II VI we identify the different 
regimes of robustness of AQC against decoherence, and 
we analyze the optimal operation conditions for a proto- 
typical AQC algorithm in the presence of decoherence. In 
Sec. |V]we discuss the scaling of decoherence effects with 
system size. Sec. I VII presents some concluding remarks. 



II. LANDAU-ZENER PROBLEM WITHOUT 
DECOHERENCE 



We start our discussion by introducing a prototypi- 
cal example of an AQC algorithm, namely the basic LZ 
problem. We therefore consider a two state system, and 
we use the spin 1/2 language, where the two states are 
called |t) and ||). In the absence of coupling to the en- 
vironment, we take the time-dependent Hamiltonian: 



A „ vt ^ 



(1) 



where A/2 is the tunnelling matrix element between the 
states |t) and ||), v is the sweep rate of the energy bias 
between the two states, and a a are the Pauli spin matri- 
ces. The instantaneous (i.e., adiabatic) two- level energy 
spectrum as a function of vt is schematically shown in the 
inset of Fig. 1. Note that the ground state and the ex- 
cited state at the degeneracy point (given by vt = 0) are, 
with the proper phase definitions, the symmetric and an- 
tisymmetric superpositions of the eigenstates evaluated 
very far from the degeneracy point. If the system is ini- 
tially in its ground state at t — > — oo, the probability for 
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the system to end up in the new ground state at t — > oo 
is given by 

P L z = l-exp{-^}. (2) 

In particular, if the system crosses the degeneracy region 
extremely slowly (v — > 0), the system is guaranteed to 
end up in the new ground state. From now on, we shall 
refer to the probability that the system ends up in the 
new ground state as the success probability, since that 
situation represents a successful run of this prototypical 
AQC algorithm. 



III. SCALING OF MINIMUM GAP WITH 
SYSTEM SIZE 

Before going into any details regarding decoherence, it 
is worth mentioning here one of the most relevant open 
questions in the study of AQC, namely the dependence 
of the minimum gap between the ground state and first- 
excited state on the system size [lOj . Since the size of 
that gap sets an upper bound on the allowed sweep rate, 
an increasingly small gap could deem an AQC algorithm 
ineffective to solve a given problem, especially in the case 
of an exponentially decreasing gap. Although that sce- 
nario would also make the algorithm more susceptible to 
decoherence, the scaling of the gap is not directly related 
to the present discussion. We shall therefore not dwell 
upon that question in this paper, and we shall leave any 
dependence of the minimum gap on system size implicit. 
Incorporating a given dependence into our results can be 
done straightforwardly. 



IV. LANDAU-ZENER PROBLEM WITH 
DECOHERENCE 

Let us start by presenting an argument that is some- 
times used to suggest robustness of AQC against deco- 
herence. We divide the noise effects into high-frequency 
and low-frequency contributions. High-frequency noise 
is responsible for relaxation processes (i.e., transitions 
between different energy levels), whereas low- frequency 
noise is responsible for dephasing processes. If we as- 
sume that the temperature is lower than the minimum 
gap encountered while running the algorithm [l|, the 
excitation rate will always be small in comparison to the 
de-excitation rate, and the system will relax to the new 
ground state at the end of every LZ crossing if neces- 
sary. High-frequency noise can therefore be neglected. 
Now, since the system is always in an eigenstate of the 
Hamiltonian, namely the ground state, dephasing is irrel- 
evant. Low- frequency noise, which describes dephasing 
processes, can therefore be neglected as well. One would 
therefore conclude that AQC is robust against decoher- 
ence. 



Given that the above argument gives strong support 
to AQC over SQC, we now discuss in some detail its 
applicability in different possible situations. An impor- 
tant point to note here is that the argument implicitly 
uses perturbation-theory results regarding relaxation and 
dephasing processes. That approach is valid only when 
the noise amplitude is small compared to the qubit en- 
ergy scales. In particular, if the assumption of small 
amplitudes in the noise signal is abandoned, the argu- 
ment breaks down. As we shall discuss in Sec. [V] this 
breakdown seems to be the case for a scalable AQC sys- 
tem. Furthermore, relaxation between macroscopically 
distinct quantum states after the LZ crossing should be 
negligible. 

A number of different approaches have been used to 
study the effects of decoherence on the LZ transition 
probability [H, 0, [H, El, E3 • Although those approaches 
are based on different underlying assumptions, they all 
produce similar qualitative results (note that they have 
different predictions regarding certain details). In partic- 
ular, all of them predict the possibility of having a max- 
imum in the success probability as a function of sweep 
rate (see Fig. 1). 

Since we shall treat a number of qualitatively differ- 
ent cases, it would be difficult to use a single model to 
describe the effects of the environment on the success 
probability. We shall therefore use two different models: 
one with a classical noise signal and one with an envi- 
ronment of quantum modes. In addition, we shall use 
thermodynamics principles when necessary. 

Before analyzing the effects of the environment on the 
system dynamics, we must specify the system operator 
involved in the system-environment coupling. In the sim- 
ple two-level problem that we are considering, that oper- 
ator must be one, or a combination, of the Pauli matrices, 
assuming the coupling is described by a product of a sys- 
tem operator and an environment operator. We note that 
away from the crossing region coupling through the oper- 
ator a z only causes dephasing, whereas coupling through 
the operators a x and a y causes relaxation. In a macro- 
scopic system, relaxation processes between macroscopic 
states are generally exponentially small. We therefore 
approach the problem at hand by taking the system in- 
troduced in Sec. [Til an( l adding a decoherence term that 
couples to the system through the operator a z . Although 
in general more complex models (i.e., many-level models) 
must be used to obtain a more detailed description of the 
effects of noise on a large AQC system, the arguments 
given below provide an initial understanding of some of 
the main mechanisms involved in the problem. 



A. Classifying the noise according to amplitude 
and correlation time 

We divide our discussion into four cases, determined 
by the following procedure: we take a noise signal with 
characteristic amplitude A (in energy units) and corre- 
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FIG. 1: Success probability P, i.e. the probability to end up in 
the new ground state after a Landau-Zener crossing, as a func- 
tion of A 2 /«, where A is twice the tunnelling matrix element, 
and v is the energy-bias sweep rate. The solid line corresponds 
to the case of no decoherence. The dashed lines correspond 
to the case of intermediate levels of decoherence (essentially 
using the classical-noise model) ; the curves were obtained fol- 
lowing Ref. [TJJ with dephasing rate T 2 (t — » ±oo) = A/200, 
A/20, and A/5. The dotted line corresponds to the limit of 
infinitely strong decoherence. Inset: Schematic view of the 
instantaneous two-level energy spectrum as a function of the 
energy bias vt. 



lation time t. We note that the power spectrum of the 
noise signal would be characterized by a (maximum) fre- 
quency Wmax that is related to the correlation time r by 
w max = 1/t. The noise power spectrum is then of order 
A 2 /uj max up to frequency w max and decreases to zero at 
higher frequencies. Note also that if the noise signal has a 
non-zero average value, we define A as the deviation from 
that average value. Depending on whether A is smaller or 
larger than the gap A, the noise is characterized as low- 
or high-amplitude noise. Similarly, depending on the re- 
lation between w max and A, the noise is characterized as 
having short or long correlation time. 



1. Low- amplitude noise with short correlation time 

We start with this case because it allows the use of 
the simple perturbation-theory results mentioned above. 
We focus on relaxation processes, because pure dephas- 
ing processes cannot have a larger effect than relaxation 
processes (note that relaxation dynamics automatically 
contains dephasing), and therefore including those can- 
not change the qualitative results we shall give below. We 
also neglect de-excitation processes for a moment. Away 
from the degeneracy region, the transition rate from the 
ground state to the excited state is negligible because 
the noise couples to the system through the operator a z , 



which is almost parallel to the system Hamiltonian. We 
therefore focus on the dynamics when the system is close 
to the degeneracy point. Since the noise power spectrum 
extends to frequencies higher than A, one finds the ex- 
citation rate from the ground state to the excited state 
around the degeneracy point to be 



r - 



A 2 



(3) 



which is essentially the noise power spectrum at the tran- 
sition frequency. One therefore straightforwardly finds 
that the time spent traversing the LZ crossing must be 
shorter than l/ro-»i ~ oJm&x/A 2 if the noise effects are 
to be minimized. Combined with the condition that the 
traversal time must be larger than 1/A, one can deter- 
mine the ideal range of sweep rates for optimal AQC 
operation. If we take the noise-induced excitation prob- 
ability to be 



A 2 A 



excited by noise 



and the LZ transition probability to be 

f ^ A2 

"excited by LZ ~ exp < ^— 



(4) 



(5) 



and we minimize the sum of those two terms, we find 
that the optimal value of v is roughly given by 



A 2 



^optimal 



ln(w max A /A 2 ) ■ 



(6) 



Similarly, one can estimate that the maximum achievable 
success probability will be 1 — Pf a ii U rc, with 



A 2 



failu 



(7) 



Note that the optimal sweep rate ^optimal depends loga- 
rithmically on the noise amplitude. That result implies 
that ^optimal can be only a few times smaller than A 2 
even if the noise power spectrum is orders of magnitude 
smaller than A. 

We note here that if one is considering the case where 
the temperature ksT is smaller than the gap A, the ex- 
citation rate will be smaller than the de-excitation rate 
by a factor exp{— A/fe^T}, and the thermal-equilibrium 
occupation probability of the excited state is given by 
1/(1 + exp{A//csT}). Therefore the above results ap- 
ply only if the expression in Eq. is smaller than 
the thermal-equilibrium occupation probability. Oth- 
erwise, one must take the de-excitation rate into ac- 
count. One then finds that the maximum obtainable 
success probability is given by the thermal-equilibrium 
value 1/(1 + exp{— A/fc^T}), and it is achieved using a 
slow sweep such that thermal equilibrium is reached. 
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2. Low- amplitude noise with long correlation time 

Since the noise amplitude is small, one can still think of 
the noise effects in terms of the transition rate from the 
ground state to the excited state. The transition rate 
in this case can be thought of as a high-order process 
[l3| . Thinking of the noise as a harmonic-oscillator bath, 
we find that an nth-order process is required to excite 
the two-level system, with n — Int(A/cj max ) + 1, and 
the function Int(x) gives the highest integer smaller than 
x. For a more concrete visualization, one can think of a 
photon bath, such that the sum of n photon energies is 
required to excite the system from the ground state to 
the excited state. The transition rate is therefore 



rv 



.4 2 



(8) 



The above expression for the excitation rate suggests that 
for the noise-driven excitation probability to be negligi- 
ble the time taken to traverse the LZ crossing must be 
smaller than l/r _>i ~ (w max /A 2 ) x (A/A) 1 " 2 ™. Given 
that A is smaller than A, the upper bound on crossing 
time above is much larger than 1/A. This case is there- 
fore the ideal case for performing AQC, allowing a high 
success probability when an appropriately chosen sweep 
rate is used. An estimate of the optimal sweep rate and 
the maximum achievable success probability can be ob- 
tained similarly to what was done in Sec. IIV A 11 In this 
case one finds 



^optimal 
urc 



A 2 



ln(w max 
A 2 

vvn ax A 



A 2n/^2n+l) 

A- 2n -' 
A 



(9) 
(10) 



Note that the characteristic noise frequency w max can- 
not be larger than the temperature fcsT, so that the 
lowest possible value of n is roughly 



Int 



A 



(ii) 



Note also that if the expression for Pf a n U re above is larger 
than 1/(1 + cxp{A/fcsT}) the optimal approach would 
be a slow sweep such that thermal equilibrium is reached 
during the crossing. 



3. High-amplitude noise with long correlation time 

We now take a slowly varying classical noise signal with 
an amplitude larger than A (note that the slowness is de- 
termined by comparison to the inverse of the gap) . We 
also take the system to be biased close to or at the de- 
generacy point. Since the amplitude of the noise signal is 
larger than the gap, one cannot use perturbation-theory 
results to describe transitions between the different eigen- 
states. Instead, one can now think of the noise signal as 



repeatedly driving LZ crossings, with noise-driven sweep 
rate 



.4 



A UJr, 



(12) 



The LZ transition probability [1 — Plz, with Plz given by 
Eq. $2$] with sweep rate v cnv is therefore not necessarily 
small, even if w max is much smaller than the gap. In 
particular, the transition probability in an environment- 
induced LZ crossing is (very roughly) given by 



P 



excited by env.— ind. LZ 



exp 



7rA 2 T 

2A 



(13) 



Given enough time, the system will therefore reach a 
state where both eigenstates have equal occupation prob- 
abilities. However, because of the exponential depen- 
dence of the transition probability on the noise param- 
eter, one can say that if the condition itA 2 t/2A ^> 1 
is satisfied, the environment-induced LZ transition prob- 
ability will be small enough that a high success rate is 
always achievable with a properly chosen value of v. The 
above criterion therefore provides the condition for high- 
amplitude noise to have a negligible effect on the success 
probability. 

One might now raise the following possibility: taking a 
LZ situation where the parameters are swept across the 
degeneracy region, one can estimate that the number of 
noise-driven crossings is of the order of A/vt. There- 
fore, if the sweep rate v is substantially larger than A/t, 
no environment-driven crossings will occur, suggesting 
that it might be possible to avoid environment-driven 
LZ transitions even if the condition wA 2 t/2A 3> 1 is not 
satisfied. It is straightforward to verify, however, that in 
order to do so one would require a value of v larger than 
A 2 . That situation would result in a high bias-driven 
LZ transition probability and, therefore, a low success 
probability. 



4- High-amplitude noise with short correlation time 

In this case one can follow the above arguments for the 
high-amplitude, low-frequency noise. Using the expres- 
sions of Sec. IIV A3! one immediately finds the intuitively 
obvious result that the success probability is 50% for low 
sweep rates and is smaller than that value for fast sweep 
rates (see the dotted line in Fig. 1). Note that the value 
50% describes the case where the two eigenstates have 
equal occupation probabilities at the end of the process. 
Note also that since we have in mind macroscopic states, 
we neglect the possibility that the system could relax to 
the ground state long after the LZ crossing. 



V. SCALABLE SYSTEM 

We now turn to the question of how decoherence effects 
scale with system size in an AQC setting with a large 
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number N of qubits (we use the typical picture of two- 
state qubits). 

We have discussed in Sec. IIVI that for large-amplitude 
noise one must think of different decoherence mechanisms 
than the usual perturbation-theory relaxation and de- 
phasing mechanisms. We therefore consider the question 
of how the noise amplitude scales with system size [ijj . 
In relation to that discussion, it is useful to classify LZ 
crossings according to the number of qubits that change 
their state during the transition. That criterion is re- 
lated to, but clearly distinct from, the question of quan- 
tifying how macroscopic a quantum state is. There has 
not been any unambiguous and universally accepted for- 
mulation of such a quantity. Following Ref. [18| we use 
a common-sense definition rather than trying to formu- 
late an operational one, which seems to be a formidable 
task. The definition is then relatively simple: a given 
LZ crossing can be referred to as an M-qubit crossing 
if M qubits change their state with the other qubits in 
the system experiencing negligible changes. We can then 
speak of few-qubit and many-qubit crossings. The former 
refers to LZ crossings of the A^-qubit system where only 
a few (say, up to four) qubits change their state, even if 
the total number of qubits in the system is macroscopic. 
The other type of LZ crossings that can occur during the 
operation of an algorithm are many-qubit crossings. In 
those crossings the number of qubits that change their 
state is of order N. 

In order to demonstrate the above-mentioned distinc- 
tion between classifying quantum states and classifying 
LZ crossings, take the plausible scenario of AQC where 
one starts with a quantum state that contains negli- 
gible multi-qubit entanglement and reaches a quantum 
superposition of macroscopically distinct states during 
the calculation. Although the quantum state becomes a 
macroscopic one, it is not necessarily the case that any 
many-qubit crossings must have been encountered (think 
for example of a macroscopic quantum state generated 
by repeatedly performing two-qubit CNOT gates). One 
should also note that even if the system is in a superposi- 
tion of macroscopically distinct states, it can still undergo 
few-qubit LZ transitions. Those transitions would most 
likely occur in one or some of the branches corresponding 
to the different macroscopic states. 

We now take an M-qubit LZ crossing. The size of the 
degeneracy region is of the order of the gap A. In a 
system with a large number of degrees of freedom, one 
can still say that the crossing region is defined by being 
within distance (in units of bias parameters) A in the 
relevant M directions from the degeneracy point, i.e. the 
point where the gap takes its smallest value along the 
path of the AQC algorithm. 

If the noise signal on a single qubit moves the system 
away from the bias point by a distance of order 5, the 
sum of the noise signals acting on the M qubits moves the 
system away from the bias point by a distance of order 
\fM5. We now take a system at or near the degeneracy 
point. If the total deviation caused by the noise is smaller 



than the width of the crossing region, which is of the 
order of A, we can use the arguments of Sec. [IV] to say 
that low frequency noise can be neglected in the sense 
that it cannot excite the system from its ground state. 

In the opposite case, i.e. when the amplitude of the 
total noise signal is larger than A, one has to worry about 
environment-driven LZ transitions. Using the results of 
Sec. IIVI we find that a rough estimate of the probability 
for the system to be excited from its ground state during 
a single typical (environment-driven) crossing is given by 



P' excited by env.— ind. LZ ^ exp 



7rA 2 T 



2VMS 



(14) 



Note that the exponential dependence of the excitation 
probability on the noise signal means that the above ex- 
pression should be thought of as an optimistic estimate; 
the true excitation probability will probably be higher, 
depending on the temporal behaviour of the noise signal. 
Using the results of Sec. IIVI the criterion on the tolerable 
single-qubit noise can now be given by 



5< 



'M 



(15) 



The probability that the noise signal will excite the 
system out of its ground state therefore depends on the 
typical value of M characterizing the LZ crossings that 
are encountered during the algorithm. Given the scaling 
of the excitation probability with M, it is highly desir- 
able to follow a path in the many-dimensional parame- 
ter space such that many-qubit LZ crossings are avoided. 
This principle can therefore remain as a major considera- 
tion in designing AQC algorithms, even if the minimum- 
gap problem discussed in Sec. IIIII is solved. 

It is not clear whether in a general AQC problem a 
path that avoids all many-qubit LZ crossings exists. The 
3-SAT problem, which is a commonly studied potential 
application of AQC [19| , provides an example where it 
seems impossible to find such a path. In that problem 
one looks for a classical state of the qubits such that a 
large number of 3-qubit logical conditions are satisfied, 
e.g. the Boolean condition "(qubit 5 and qubit 24) or 
qubit 57" . In the plausible scenario where one configura- 
tion satisfies all the logical conditions but a large number 
of other, macroscopically distinct configurations violate 
only a few conditions, a quantum superposition involving 
a large number of macroscopically distinct configurations 
must be retained until near the end of the calculation, as 
they are eliminated slowly with the testing of more and 
more conditions. At that point it would require a many- 
qubit LZ crossing to eliminate those last surviving near- 
solutions in favor of the unique solution of the problem. 
The 3-SAT problem therefore appears to be one where 
decoherence can be a major obstacle. The fact that the 
path of an AQC algorithm is designed without knowing 
the quantum state that will exist at each point in the 
algorithm raises similar doubts about the possibility of a 
priori guessing the best path to follow in a general prob- 
lem. 
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The above arguments therefore raise questions that 
must be answered in designing an AQC approach in a 
macroscopic setup. Until those questions are answered, 
it is not clear to what extent AQC is less susceptible to 
noise than SQC, especially given the condition that we 
found above requiring the noise signal to decrease with 
increasing system size [l2j . 

Even if achieving the ground state is not possible, e.g. 
because of decoherence or a small minimum gap, a new 
promising proposal notes that finding a near-solution 
can, under certain conditions, be considered a success 
of the algorithm [20j • Because a high success probability 
(in the sense of Sec. [TTJ) is not required, that approach 
could be more robust against decoherence. 

It is also worth noting here that we have used the sim- 
ple model of a two-state LZ problem, which represents 
a prototypical AQC algorithm. The number of degrees 
of freedom in a large AQC setup increases with system 
size. More complex models will be required in order to 
both analyze the effects of noise and determine the op- 
timal path in those many-dimensional problems. Reach- 
ing a better understanding of the structure of the energy 
manifolds in these many-dimensional systems is therefore 
highly desirable. 

VI. CONCLUSION 

We have analyzed the effects of noise on a prototyp- 
ical AQC algorithm, namely the LZ problem. We have 



found general principles that determine the robustness 
of the algorithm against noise sources with a variety of 
properties according to their amplitude and correlation 
times. We have also determined the ideal operation con- 
ditions that are required to maximize the success prob- 
ability, and we have analyzed the scaling of noise effects 
with system size. Our results provide guidelines for the 
optimal implementation of an AQC algorithm and raise 
questions that must be answered before determining the 
suitability of AQC to tackle a given problem. Given the 
promise of AQC as an alternative approach to achieve 
extremely high-speed computation, we believe that our 
results will contribute to a better understanding of that 
approach, towards which initial experimental steps have 
already been taken [2l|, [22j . 
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